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We examine the role of thermal fluctuations in binary condensate mixtures of dilute atomic gases. In partic¬ 
ular, we use Hartree-Fock-Bogoliubov with Popov approximation to probe the impact of non-condensate atoms 
to the phenomenon of phase-separation in two-component Bose-Einstein condensates. We demonstrate that, in 
comparison to T = 0, there is a suppression in the phase-separation of the binary condensates at T 7 ^ 0. This 
arises from the interaction of the condensate atoms with the thermal cloud. We also show that, when T yl 0 
it is possible to distinguish the phase-separated case from miscible from the trends in the correlation function. 
However, this is not the case at T = 0. 


PACS numbers: 03.75.Mn,03.75.Hh,67.60.Bc,67.85.Bc 

Phase-separation in a two-component fluid is ubiquitous in 
nature, and the transition from miscible to immiscible phase 
is a quintessential example of critical phenomena. One clas¬ 
sic example is the temperature driven phase-separation in the 
cyclohexane-aniline mixture [1]. It is then natural to ask what 
are the similarities and differences in binary mixtures of quan¬ 
tum fluids ? The recent experimental advances in binary Bose- 
Einstein condensates (BECs) of dilute atomic gases provide 
an ideal testbed to address such a question. In the case of bi¬ 
nary mixtures of BECs or two-species BEC (TBEC) tuning 
the interaction through Feshbach resonances [2, 3] can render 
it miscible or immiscible. Using improved experimental tech¬ 
niques, TBECs have been achieved in mixtures of two differ¬ 
ent alkali atoms [4-8], or two different isotopes [2, 9, 10] and 
atoms of the same element in different hyperfine states [ 1 1 - 
14] over the last decade. The remarkable feature of phase- 
separation in TBECs has been successfully observed in S 5 Rb- 
s, Rb [9, 10] and 87 Rb- 133 Cs [7] condensate mixtures. 

The criterion for phase-separation, derived from Thomas- 
Fermi (TF) approximation at zero temperature [15], is that the 
intra-(f7n, C/ 22 ) and inter-species interaction (C/ 12 ) strengths, 
must satisfy the inequality Uf 2 > C /11 C/ 22 . However, experi¬ 
ments are conducted at finite temperatures, and therefore, de¬ 
viations from the criterion is to be expected. Theoretical stud¬ 
ies on effects of thermal cloud on phase-separation have been 
carried out for homogeneous binary Bose gases using Hartree- 
Fock theory [16] and Iarge -N approximation [17], Phase- 
separation of trapped binary mixtures at finite temperature has 
also been examined using local-density approximation [18]. 
In this Letter we address this issue by using Hartree-Fock- 
Bogoliubov theory with Popov approximation (HFB-Popov) 
[19] to account for the thermal fluctuations. It is a gapless 
formalism satisfying Hugenholtz-Pines theorem [20] and can 
be employed to compute the energy eigenspectra of the quasi¬ 
particle excitations of the condensates. 

The method has been validated extensively in single species 
BEC, and we have used it in our recent works to examine the 
effect of quantum fluctuations in TBECs [21]. In the present 
work, we systematically study the role of thermal fluctuations 
in the phenomenon of phase-separation in trapped TBECs. 
Our studies reveal that at T / 0, the constituent species in the 
TBEC undergo phase-separation at a higher U 12 than the value 
predicted based on the TF-approximation at T = 0. Consis¬ 


tent with experimental observations of dual species conden¬ 
sate of s, Rb and 133 Cs [7], our theoretical investigations show 
that even when the phase-separation condition is met, there is 
a sizable overlap between the two species. We attribute this 
to the presence of the thermal cloud, which have profound af¬ 
fects on the miscibility-immiscibility transition. At T = 0, 
the TBECs are coherent throughout the spatial extent of the 
condensate, however, when T 7 / 0 coherence decays and is re¬ 
flected in the correlation function. This implies that at T = 0, 
the miscible or immiscible phases are indistinguishable from 
the trends in the correlation function. But, for T 7 ^ 0 the 
miscible-immiscible transition and the associated changes in 
the density profiles have a characteristic signature in the form 
of the correlation functions. There is a smooth cross-over 
between the correlations functions when the transition oc¬ 
curs. Interspecies Feshbach resonances of ultracold bosons 
have been experimentally demonstrated for Na-Rb [22], K- 
Rb [23] and Cs-Rb [24] mixtures, but Bose condensed mix¬ 
ture of Na-Rb is yet to observed experimentally. The Cs-Rb 
condensate mixture is a stepping stone towards the produc¬ 
tion of quantum gas of dipolar RbCs molecules, as unlike 
the KRb molecule, the rovibrational ground state of RbCs 
molecule is stable against exchange of atoms. Considering 
this we focus our study on the finite temperature effects in 
the Cs-Rb condensate mixture. Other than tuning the inter¬ 
species, it is also possible to steer the condensate mixture 
through the miscible-immiscible transition using intraspecies 
Feshbach resonance. An example is the tuning of the in¬ 
traspecies interaction of 85 Rb in 85 Rb- 8 'Rb [9, 10], and for 
this system too, we have examined the suppression of phase- 
separation at T 7 ^ 0 [25], It must be emphasized that, as the 
background scattering length of So Rb is negative, it is possi¬ 
ble to obtain 85 Rb BEC [26] only with the use of Feshbach 
resonance [27]. 

Theory — We consider a cigar shaped TBEC, where the 
frequencies of the harmonic trapping potential satisfy the con¬ 
dition uj± u z with ui x = u} y = oj±. In this case, the 
radial excitation energies are large and assume the radial de¬ 
grees of freedom are frozen for which hoj± Pk- So, the 
dynamics and hence the excitations occur only along the axial 
direction, z-axis, of the trap. In the mean-field regime, using 
HFB-Popov approximation [21, 28], a pair of coupled gen¬ 
eralized ID Gross-Pitaevskii (GP) equations describe the dy- 
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namics and density distribution of the TBEC. The combined 
form of the equations is 

h k <j) k + u kk [n ck + 2 h k \ <p k + Ui 2 n^~ k (j) k = 0, (1) 

where h k = (— h 2 /2m k )d 2 /dz 2 + V k (z) — n k is the one-body 
part of the Hamiltonian, with k = 1,2 as the species label. 
The strength of the coupling constants are given by U kk = 
{o, kk \)/m k and t / 12 = (ai 2 A)/( 2 toi 2 ), where, for cigar 
shaped traps A = ui±/uj z 1. Without loss of generality, for 
stable configurations the intra-species scattering lengths a kk , 
and the inter-species scattering length ai 2 are considered as 
positive (repulsive). Under the HFB approximation, the Bose 
field operators are decomposed as T/,. = <p k + yj k , where the 
cf> k s are the stationary solutions of Eq. (1) obtained by evolv¬ 
ing the solution in imaginary time, with n ck (z) = \<p k (z)\ 2 . 
The field operator ip k {z) represents the fluctuation part of 
( 2 ), it incorporates both the quantum and thermal fluctua¬ 
tions. The fluctuation operator, both quantum and thermal, are 
functions of the elementary excitations of the system, which 
solves the coupled Bogoliubov-de-Gennes equations 

Cyuij - Unplvij + U 12 P 1 (4>2 u 2j - 4>2V2j ) = EjUij, (2a) 

LiVij + Unp^uij - {<t>2V2j - <p\u 2 j) = EjVij, (2b) 

t 2 u 2 j - U 2 2 (fiv 2 j + U 12 P 2 - (filVij) = EjU2j, (2c) 

£ 2 V2j + U 2 2(pl 2 U2j - U 12 PI {4>lVij - = EjV 2 j, (2d) 

where C\ = (h\ + 2Unn\ + f/i 2 n 2 ), C 2 = (/i 2 + 2f7 22 ro 2 + 
U^n-i ) and C k = —C k . Here u k jS and v k j s are the Bogoli- 
ubov quasi-particle amplitudes corresponding to the jth en¬ 
ergy eigenvalue. The quantities h k (z) = ( ■ip\(z,t)'ip k (z,t )), 
and n k (z) = n ck (z) + h k (z) are defined as non-condensate, 
and total density, respectively. To solve the above eigenvalue 
equations, the u k j s and v k) s are decomposed into a linear 
combination of harmonic oscillator eigenstates. The order pa¬ 
rameters (j) k s and the non-condensate densities fi k s are then 
the self-consistent solutions of the coupled Eqns. (1) and (2). 
The thermal components, in terms of the quasi-particle ampli¬ 
tudes, is 


n k = EUK'I 2 + M 2 )MEj) + \v kj \ 2 }, (3) 
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where, N 0 (Ej) = (e^ Ej — 1 ) _1 with /3 = 1 /(/cbT) is the 

Bose factor of the jth quasi-particle mode at temperature T. 
A more detailed description of the decomposition and deriva¬ 
tion of the relevant equations are given else where [25]. In this 
Letter we examine the role of temperature in phase-separation 
of TBECs. For this, a measure of phase separation is the over¬ 
lap integral 

= [/ni(z)n 2 (z)rfz ] 2 

[f n i(z)dz] [jnl{z)dz ]' 

Miscible phase is when A = 1 and signifies complete over¬ 
lap of the two species, whereas the binary condensate is com¬ 
pletely phase-separated when A = 0 [29], 


T = 0 T = 5 T=10 



FIG. 1. (Color online) The suppression of phase-separation in 
87 Rb- 133 Cs TBEC at a 12 = 295a 0 . (a)-(c) The solid and dashed 
red (black) lines represent ncs(ftRb) and ncs(uRb), respectively, 
at T = 0, 5, lOnK. (d)-(f) The solid red (black) lines represent 
n c Cs(rtcRb) at T = 0, 5, lOnK respectively. The dashed red (black) 
lines n c Cs(n C Rb) at T = 0 with the same number of condensate 
atoms at T = 0, 5, lOnK respectively. Here, n and 2 are measured 
in units of a]7 a c and a oac , respectively. 


In terms of the Bose field operator d/ k , the normalized first 
order or the off-diagonal correlation function, which is also a 
measure of the phase fluctuations, is 

(^(2)^(2')) (5) 

(^(2)T- fc (2))(T'i(20^ fc (^)} 

It can also be expressed in terms of off-diagonal condensate 
and non-condensate densities as 

( 1 )/ /x n ck (z,z') +h k (z,z') 

9 k (z,z) = - ; (6) 

V n k {z)n k \z') 



where, 

n ck {z,z') = 4>* k {z)(j) k {z') 

n k (z,z') = Y2{[u k j(z)ukj(z') +vl j {z)v kj (z')}N 0 {E j ) 
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+v* k j{z)v kj {z')}. 

At T = 0, when the entire system is coherent and character¬ 
ized by the presence of condensate only, then = 1 within 
the extent of the condensate, whether it is in the miscible or 
in the immiscible regime. So, one cannot distinguish between 
the two phases from the nature of the correlation functions of 
the individual species. However, at T ^ 0, a clear signature 
of miscible-immiscible transition of the density profiles is re¬ 
flected in the form of the correlation functions. 

The thermal suppression of phase-separation is generic to 
any binary condensate mixture. However, for comparison 
with experimental realizations we consider the 133 Cs- 8l Rb 
BEC mixture with 133 Cs labelled as species 1, and Sl Rb as 
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species 2. Here after, for brevity we drop the mass numbers 
and write these as Cs and Rb. The intra-species scattering 
lengths are an = ac s = 280a o , a 22 = aRb = 100 a o , the 
inter-species scattering length is ai 2 = acsRb = 295ao with 
Ncs = -^Rb = 5 x 10 3 , and ao as the Bohr radius. To form 
a quasi-ID trap we take u> z (c s ) = 27r x 4.55Hz and w z (R b ) = 
27 t x 3.89Hz; Wj_(c s ) = 50w 2 (c s ) and Wj_( Rb ) = 50w 2 ( Rb ). 
For this value of ujj_, the temperature along the radial direc¬ 
tion is hoj±/kB ~ llnK, and the tight confinement condition 
is valid as gk/tujJ± ~ 10~ 2 . In addition to this, the heal¬ 
ing length 1 /n/ :: . Thus the system is in the weakly 

interacting TF regime [30] and mean field description through 
GP-equation is valid. For this parameter set, the ground state 
density distribution is phase-separated with species 1 at the 
center and surrounded by the species 2 at the edges. We re¬ 
fer to this configuration of density profiles as sandwich type. 
This choice of parameters is consistent with the experimen¬ 
tal parameters of a recent work on quasi-ID TBEC of differ¬ 
ent hyperfine states of tS 'Rb [31], in which dynamical evo¬ 
lution of mixtures of quantum gases has been observed. It 
should be emphasized here that sandwich type density pro¬ 
files are applicable only to trapped systems. In uniform sys¬ 
tems, at phase-separation, the energetically preferred states 
are the symmetry-broken density profiles where one species 
is entirely to the left and the other is entirely to the right. We 
refer to this configuration of density profiles as side-by-side 
type. In the present work, we demonstrate the role of ther¬ 
mal cloud in sandwich type density profiles since these are 
unique to trapped systems and experimentally pertinent. For 
the homogeneous binary condensates, using periodic bound¬ 
ary condition with w 2 = 0 in our computations, we do get 
side-by-side density profiles at phase separation. As an exam¬ 
ple, the density profiles are shown in Ref. [25] and these are 
consistent with the results reported in previous works [32], In 
the computations the spatial and temporal variables are scaled 
as 2 /a osc (cs) and u> z (c s )t to render the equations dimension¬ 
less. 

At T = 0, in TBECs, as mentioned earlier, the criterion 
for phase separation is U12 > \JV\ i U22. With the param¬ 
eters of Cs-Rb TBEC, consider keeping ac s and aR b fixed, 
but varying ai 2 = «CsRb through a magnetic Feshbach res¬ 
onance [24], The condition for phase-separation, using TF- 
approximation, is then ai 2 > 261ao. When ai 2 = 0, the 
TBEC is non-interacting and the two species are completely 
miscible, in which case A = 1. On increasing ai 2 , the ex¬ 
tent of overlap between the two-species decreases, and hence 
A decreases. For instance, at ai 2 = 50ao, A = 0.97 and 
it decreases monotonically with A —X 0 at complete phase- 
separation. At ai 2 = 295ao, just at the onset of phase- 
separation, A = 0.14. As shown in Fig. 1(a), the density pro¬ 
files corresponding to the two-species have interfacial overlap, 
and the interaction parameters satisfy the phase-separation 
condition. Furthermore, at phase separation, n c c s (0) is maxi¬ 
mum, whereas n C R b (0) ~ 0 and the species do not have sig¬ 
nificant overlap. In other words, Cs at the center of the trap 
is flanked by Rb at the edges and A ss 10 _1 . It is also to be 
mentioned here that for phase-separation, there is consider¬ 
able difference between the values of ai 2 derived from TF ap¬ 


proximation and the numerical solution of GP equation. This 
can be attributed to large gradients in condensate densities, 
which are ignored in the TF-approximation. 


T = 0 T = 5 T= 10 



FIG. 2. (Color online) (a)-(c) The first order spatial correlation func¬ 
tion, 3cs/Rb(,0> z ) with z ^ 0, of 8 ' Rb- 133 Cs TBEC at equilibrium 
for ai 2 = 295ao at T = 0, 5, lOnK respectively. Here 2 is measured 
in units of a osc . 

Suppression of phase segregation — For T / 0 the Bose 
factor No f 0 , so in addition to the quantum fluctuations, 
the non-condensate densities have contributions from the 
thermal cloud as well. The condensate atoms n c k then in¬ 
teract with fik of both the species, and modify n ( ± of both the 
species. For illustration, at T = 5nK and « 12 = 295ao, the to¬ 
tal and non-condensate density profiles are shown in Fig. 1(b). 
Compared to the density profiles in Fig. 1(a), there is a re¬ 
markable change in n c Rb as a result of the finite temperature: 
ficRb(O) > 0. Thus, keeping all the parameters same, but 
taking T = 5nK, the two species have substantial overlap as 
shown in Fig. 1(b) and A becomes « 0.55. In other words, 
the finite temperature transforms the phase separated TBEC at 
T = 0 to a partially miscible phase. The degree of overlap in¬ 
creases with temperature and at T = lOnK, the TBEC is mis¬ 
cible as A « 0.77. Thus, with the increase in temperature, the 
density of thermal cloud increases and the phase-separation 
is suppressed. This is evident from Fig. 1(c), which shows the 
plots of corresponding total and non-condensate density pro¬ 
files. Thus, ai 2 has to be greater than 295ao at T ^ 0 for 
phase-separation to occur. To confirm that the suppression is 
a consequence of non-zero temperature, we identify and com¬ 
pute the number of condensate atoms in each species, and use 
these numbers for T = 0 computations. Despite the differ¬ 
ence in the numbers of atoms, from the plots in Figs, l(d-f), 
the TBEC retains the immiscible profiles at zero temperature. 
This implies that without the thermal cloud, there are no de¬ 
viations from the usual phase-separation condition. For com¬ 
parison, plots of n c fc from the finite temperature cases are also 
shown in the figure. 

To investigate the spatial coherence at equilibrium, we 
examine the nature of the first order correlation function 
gP(z,z') as defined in Eq. ( 6 ). As to be expected, profile 
of (z,z r ) depends on the interplay between the interaction 
strength and temperature. As stated earlier, at T = 0, there is 
perfect coherence in the Cs-Rb TBEC and fl , 0 s / Rb (O, z ) = 1 
within the spatial extent of the condensates. This is indepen¬ 
dent of whether the TBEC is in miscible or immiscible regime. 
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FIG. 3. (Color online) (a)-(d) The first order spatial correlation func¬ 
tion, ffcs/Rb(0, z ) ~ ^ 0, of s 'Rb- 153 Cs TBEC at equilibrium 
at T = 5nK for a ±2 = 0, 220, 250, 290ao, respectively. Here z is 
measured in units of a osc - 


acsRb = 650ao 
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FIG. 4. (Color online) Density profiles showing complete phase- 
separation at T = 0 and T = lOnK. (a) Phase-separation in 8 'Rb- 
138 Cs TBEC for ai 2 = 650ao. The solid and dashed orange (brown) 
lines represent n c c s (n C Rb) and n c c s (ftcRb), respectively at T = 0. 
The solid and dashed red (black) lines represent n c cs(n C Rb) and 
KcCs (ttcRb), respectively at T = lOnK. (b) Phase-separation in 
85 Rb- 87 Rb TBEC for 012 = 20ao- The solid and dashed orange 
(brown) lines represent n c and n of 8> Rb ( s ‘Rb), respectively at 
T = 0. The solid and dashed red (black) lines represent n c and n of 
8 °Rb ( 8 ‘ Rb), respectively at T = lOnK. Here, n and z are measured 
in units of a)7sc and flosc, respectively. 


For simplicity and based on the symmetry of the system we 
consider <?c s / Rb ( 0 , z) with z ^ 0 , and plots at different tem¬ 
peratures are shown in Fig. 2. At T = 0 the form of the 
ffCs/Rb( 0 > z ) remains unchanged as the system undergoes the 
dramatic transition from miscible to immiscible phase. This 
is evident from the plot in Fig. 2(a). However, when T ^ 0, 
unlike the zero temperature case, Scs/Rb( 0 ) z ) * s maximum 
at z = 0 and decays to zero with z. This is due to the 
non-condensate atoms, which modify the the nature of coher¬ 
ence in the system. The rate of decay of the <?c s /Rb (®> z ) ^ n " 
creases with temperature, and this is evident from the plots of 
3cs/Rb(0> ~) at T = 5nK and T = lOnK shown in Figs. 2(b- 
c) for ai 2 = 295ao- We also observe a dramatic variation 
in fl , Cs/Rb( 0 ; z ) at fixed temperature, but the value of 012 is 
steered from miscible to immiscible regime. At the outset, 
when the TBEC is miscible at 012 = 0, g p b (0, z) decays to 0 
at a larger distance than g^'j (0, z) as shown in Fig. 3(a). This 
is because nR b has a larger spatial extent than nc s - As 012 is 
increased, the TBEC undergoes a phase-transition from mis¬ 
cible to sandwich type density profiles. Along with this, the 
distance at which .g bb ( 0 , z) falls off to zero increases with in¬ 
crease in ai 2 . On the contrary, the distance at which g^J (0, z) 
falls off to zero decreases with increase in 012 - This causes the 
gjj 1 ( 0 , z) of the individual species to cross each other at a cer¬ 
tain zo. At zo, the two species have equal <?c s /R b (0j zo) an d 
this is a characteristic signature of immiscible phase. These 
features are shown in Figs. 3(b-d). It deserves to be mentioned 
here that zo increases, and ffcs/Rb(*-*’ z o) decreases with in¬ 
crease in (I 12 . In addition, there is a dramatic difference in the 
decay rates of ffc s / Rb (0, Zq); it is much faster in Cs. This is 


attributed to the fact that both n cRb and n Rb increase along z 
within the bulk of Cs-Rb TBEC. Where as in Cs, around the 
origin n c c s decreases but hc s increases. This trend is similar 
with a single-species Cs condensate. The presence of Rb does 
not affect the nature of g^'j (0, z) in Cs-Rb TBEC. Around the 
point of phase-separation n C R b (0) in a Cs-Rb TBEC is dis¬ 
tinctly different from single species Rb condensate, so is the 
nature of g^ (0, z) [25]. 

Segregation independent of temperature — In the domain 
of large 012 , U 12 A rj 0, and the phase- 

segregation is more prominent. However, due to the geometry 
of the TBEC mean-field approximation is still valid. In this 
domain the interfacial overlap is minimal and the TBECs as¬ 
sume sandwich type density profile. The system is then equiv¬ 
alent to three coupled condensate fragments, and as a result, 
the Bogoliubov analysis shows the presence of three Gold- 
stone modes [28]. For the Cs-Rb TBEC considered here, the 
background inter-species scattering length ac S Rb = 650ao 
satisfies the above condition. With this value of < 212 , at T = 0 
as shown in Fig. 4(a), Cs condensate lies at the center of 
the trap and Rb condensate at the edges. So, at the center 
tiRb(O) = 0 and ncs(0) is maximum. With the increase 
in 012 , there is a decrease in the number of non-condensate 
atoms arising from quantum fluctuations. This is a manifesta¬ 
tion of smaller overlap between the condensates at the inter¬ 
faces. On the contrary, for a single species BEC, with the in¬ 
crease in intra-species interaction strength, the number of non¬ 
condensate atoms due to quantum fluctuations increases[ 21 ]. 
When 0, the thermal density hk interact with the conden¬ 
sate clouds through the intra- and inter-species interactions. 
But, due to the large < 212 , the inter-species interaction energy 
is much larger than the intra-species interaction energy. This 
makes n Rb ( 0 ) ~ 0 , and there is little overlap of the thermal 
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cloud of one species with the condensate of the other species, 
such that A < 0.1. Thus, there is no thermal suppression in 
the the domain of large a 12 . We observe similar results in the 
case of 85 Rb- 8 'Rb TBEC as well, where the the intra-species 
interaction of 85 Rb is decreased to obtain completely phase- 
separated density profiles. These are shown in Fig. 4(b). 

Conclusions — At finite temperatures, to examine the prop¬ 
erties of binary condensates in the neighbourhood of phase 
separation, it is essential to incorporate the thermal compo¬ 
nent. In general, there is a delay or suppression of phase- 
separation due to the thermal component, and we have ex¬ 
amined this in detail with the Cs-Rb binary condensate as an 
example. In this system the transition is driven by tuning the 
inter-species interaction, and similar results are obtained in 
85 Rb- 8 'Rb binary condensate, where tuning the intra-species 
interaction of 85 Rb induce the transition. The binary con¬ 
densate mixtures of dilute atomic gases are different from 
the classical binary fluids which undergo miscible-immiscible 
transition with temperature as control parameter. First, the 
variation of temperature in TBECs is applicable only below 
the lower of the two critical temperatures. Second, each 
species has two sub-components, the condensate and non¬ 
condensate atoms. The condensate or the superfluid compo¬ 
nents are coherent, but the non-condensate components are 
incoherent and like the normal gas. Third, there are spatial 


density variations of all the components due to the nature 
of the confining potential and diluteness of the atomic gas. 
Fourth, beyond a certain critical value of interaction strength 
or in the U 12 \/ U\ -\ U 22 domain, temperature does not al¬ 

ter the density profiles. Finally, the transition to the phase 
separated domain at finite temperatures is associated with a 
distinct change in the profile of the correlation function. Our 
results provide an explanation of the experimentally observed 
density profiles at phase-separation[7]. Even when the phase- 
separation condition is satisfied there is a finite overlap be¬ 
tween the two species at finite temperature which is due to the 
presence of thermal cloud. Our findings clearly demonstrate 
the dominance of thermal fluctuations over quantum fluctu¬ 
ations which causes the suppression of phase-segregation in 
any TBEC experiment. 
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